Application of a time-fractal fractional derivative with a power-law kernel to the Burke-Shaw system based on Newton's interpolation polynomials

This paper proposes some updated and improved numerical schemes based on Newton's interpolation polynomial. A Burke-Shaw system of the time-fractal fractional derivative with a power-law kernel is presented as well as some illustrative examples. To solve the model system, the fractal-fractional derivative operator is used. Under Caputo's fractal-fractional operator, fixed point theory proves Burke-Shaw's existence and uniqueness. Additionally, we have calculated the Lyapunov exponent (LE) of the proposed system. This method is illustrated with a numerical example to demonstrate the applicability and efficiency of the suggested method. To analyze this system numerically, we use a fractal- fractional numerical scheme with a power-law kernel to analyze the variable order fractal- fractional system. Furthermore, the Atangana-Seda numerical scheme, based on Newton polynomials, has been used to solve this problem. This novel method is illustrated with numerical examples. Simulated and analytical results agree. We contribute to real-world mathematical applications. Finally, we applied a numerical successive approximation method to solve the fractional model.• The purpose of this section is to define a mathematical model to study the dynamic behavior of the Burke-Shaw system.• With the Danca algorithm [1,2] and Adams-Bashforth-Moulton numerical scheme, we compute the Lyapunov exponent of the system, which is useful for studying Dissipativity.• In a generalized numerical method, we simulate the solutions of the system using the time-fractal fractional derivative of Atangana-Seda.


a b s t r a c t
This paper proposes some updated and improved numerical schemes based on Newton's interpolation polynomial.A Burke-Shaw system of the time-fractal fractional derivative with a powerlaw kernel is presented as well as some illustrative examples.To solve the model system, the fractal-fractional derivative operator is used.Under Caputo's fractal-fractional operator, fixed point theory proves Burke-Shaw's existence and uniqueness.Additionally, we have calculated the Lyapunov exponent (LE) of the proposed system.This method is illustrated with a numerical example to demonstrate the applicability and efficiency of the suggested method.To analyze this system numerically, we use a fractal-fractional numerical scheme with a power-law kernel to analyze the variable order fractal-fractional system.Furthermore, the Atangana-Seda numerical scheme, based on Newton polynomials, has been used to solve this problem.This novel method is illustrated with numerical examples.Simulated and analytical results agree.We contribute to real-world mathematical applications.Finally, we applied a numerical successive approximation method to solve the fractional model.
• The purpose of this section is to define a mathematical model to study the dynamic behavior of the Burke-Shaw system.• With the Danca algorithm [1,2] and Adams-Bashforth-Moulton numerical scheme, we compute the Lyapunov exponent of the system, which is useful for studying Dissipativity.• In a generalized numerical method, we simulate the solutions of the system using the timefractal fractional derivative of Atangana-Seda.

Specifications Table
Subject area Mathematics and Statistics More specific subject area Biomathematics Method name: Fractal-fractional Caputo Method Name and reference of original method: Atangana A, Araz SI (2021) New numerical Scheme with Newton polynomial: theory, methods, and applications.Academic, Cambridge.978-0323854481.Resource availability: This method has been developed in MATLAB

Method details
The Lorenz attractor was named after Edward Norton Lorenz, who derived it from the simplified equations of convection rolls arising in the atmosphere equations in 1963.Burke and Shaw derived the Burke-Shaw system from the Lorenz system [3] .This system has a similar algebraic structure to the Lorenz system but is topologically nonequivalent to the generalized Lorenz-type system and can be expressed as follows: ( )  = − ( x( t) + y ( t) ) , ( )  = −  x( t) z ( t) − y( t) , ( )  =  x( t) y ( t) + , where ,  , and  state variables and , , ,  are real constants.Its fractional-order version can be expressed as [4] : where  is the derivative order, which could be arbitrary real numbers.Chaos control is performed by a simple linear controller, and a numerical simulation of the control is provided.In addition, Chen, and Lee [5] introduced a novel chaotic system capable of generating dual-role chaos attractors when investigating rigid body motion anti-chaos control.Richter [6] studied the stability and chaos control of Newton-Leipnik systems [7][8][9] using static nonlinear feedback laws based on Lyapunov functions.In [9] , Long-Jye Sheu et al., investigated the dynamics of the Newton-Leipnik system with fractional order and was studied numerically.Fractional calculus (FC) allows integration and differentiation of operators in fractional order.Samko proposed a fascinating extension of constant-order FC in [10][11][12][13][14][15][16][17][18][19] .Solis-Perez, et al. [20] , introduced fractional operators that consider order as a function of time, space, or other variables.Variable-order fractional differential equations cannot be solved exactly, so developing numerical schemes for solving these equations is crucial.For fractional differential equations, the Adams-Bashforth method is highly effective -see [ 21 , 2 , 22 ].In [20] , Solis-Perez, et al. developed a constant-order numerical scheme that combines fractional calculus and Lagrange polynomials.Using this method, they generalized the numerical schemes for simulating variable-order fractional differential operators with power-law, exponentiallaw, and Mittag-Leffler kernels.See also [23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39][40] .The aim of this paper is to propose some new improved numerical schemes based on Newton's interpolation polynomial.The purpose of this study is to describe the properties of the Burke-Shaw system in the sense of a time-fractal fractional derivative with a power-law kernel.We examine the numerical aspects of the presented system as well as the existence a uniqueness of the solutions to the presented model.By using the Banach fixed point theorem and Picard iterative method, we look at the stability of fixed points and determine the range at which it can be controlled to be stable.With Lagrange polynomials, a variable-order fractal fractional with a power-law kernel is used to simulate state variables.In summary, the paper makes the following claims: (a) Demonstrate the effectiveness and uniqueness of the system presented.(b) Utilizing fractal-fractional derivatives to simulate the solution of the system using a power-law kernel method (c) Obtaining the Lyapunov exponent for the proposed system.(d) Simulating the behavior of different state variables and investigating the impact of change the values of some parameters in the model.

Model formulation
We consider the family Ψ of all increasing functions where  is a normed space.Then, ( Definition 2. [32] .Let a continuous function  ∶ (,  ) → [0 , ∞) be fractal differentiable of fractal order .Then, the fractalfractional derivative of F equipped with the power-law-type kernel of order  in the sense of Riemann-Liouville is defined by is the fractal derivative and  − 1 < ,  ≤  ∈ .Using the fractal-fractional sense of differential and integral operators we get the following Burke-Shaw model: FF  , 0 ,  (  ) =  x(t )y(t ) +  .
Since the integral is differentiable, we can rewrite the system (1) as follows Now, when we replace the derivative RL D by C D , applying fractional integral, we get the solution as follows: where
For each ,  ′ ∈ , we have This shows that F 2 is Lipschitz with respect to y with the Lipschitz constant w 2 > 0. For each ,  ′ ∈ , we have This shows that F 3 is Lipschitz with respect to z with the Lipschitz constant w 3 > 0. As a result, ,  ,  are Lipschitzian with respect to the Lipschitz constants  1 ,  2 ,  3 > 0 , respectively.

Local stability analysis
For a specific value of the parameters such as  =  =  = 10 and  = 13   = 4 .272 with the initial value (0 . 1 , 0 . 1 , 0 . 1 ) , which makes the Burke-Shaw attractor chaotic.The Jacobian matrix J of the Burke-Shaw system for the equilibrium point  * = ( * ,  * ,  * ) is defined as follows: The Burke-Shaw system (1) with above parameters has two equilibrium points: E 1 = (1 .1402 , − 1 .1402 , 0 .1) and E 2 = (− 1 .1402 , 1 .1402 , 0 .1).Obviously, due to parameter d in the model there is no equilibrium point at the origin.For the equilibrium points E 1 , E 2 , we obtain the same eigenvalues.

Kaplan-Yorke dimension
Here are the Kaplan-Yorke dimensions for some of the fractional derivatives presented in Table 1 : The fact that all the Kaplan-Yorke dimensions calculated earlier are fractional is another indication that the system is moving in a chaotic direction.Fig. 1 simulation results demonstrate the Lyapunov exponential spectrum technique for chaotic fractional-order systems high accuracy and convergence.

Dissipativity
The divergent flow of ( 1) is dissipative if and only if ∇  < 0 , If  −  < 0 .4 then the system is dissipative.
and we reorder the above equation as follows: Taking as and calculating the integral of above equation, we can get the following: and we write If we use two step Newton polynomial, we write the following: and if we organize the above equations, we can have the following.Almutairi and S. Saber MethodsX 12 (2024) 102510 We have the following calculation.
Thus, numerical solution of Burk-Shaw is given the following scheme. , Numerical scheme of variable order fractal-fractional with power-law kernel Let  be a differential function.Let  be a constant fractional order such that 0 <  ≤ 1 .Let (  ) > 0 be continuous function, then a fractional derivative of  with order  and fractal variable dimension (  ) is given by ( ) .
The new fractional integral with power-law kernel is given by In this section, we give the derivation of a numerical solution of the following problem: FF  ,( ) 0 , ( ) = h( t, u( t) ) , u(0) =  0 .Applying the new fractional integral with power kernel, we can rewrite the above equation as ( ) d .At the point   +1 = ( + 1 )Δ , we can have the following: ( ) .
Then, we have and we write Using the Lagrange polynomial, the above equation can be revised, Thus, we have And We put the calculations for the above integrals into Eq.( 4) , and we obtain the following approximation.
Thus, if we replace the function (,  (  ) ) by its value, we have the following:

Discussion
The equilibrium points of system (1) and the corresponding eigenvalues of the Jacobian matrix are shown in Table 1 .In the chaotic 3D chaos, the equilibrium points of the Burke-Shaw system (1) that yield all unstable eigenvalue as illustrated in Table 1 .A balance with exactly five unstable eigenvalues, the saddle point or saddle focus with index 2, is responsible for the generation of the rolling attractor.It is found in Table 1 that the minimum order corresponding to (1) is amin = 0.94, and chaos may exist above this order.Therefore, the theoretically calculated minimum effective size of Burke-Shaw system is 2.82, and this finding is further verified in the numerical simulation results in Section 4. The system shows better dynamic behavior.

Conclusion
Fractal-fractional operators could accurately replicate and reveal some chaos.However, because of their non-linearity, their analytical solutions are difficult to obtain and, in some circumstances, impossible to achieve due to their non-linearity.Researchers rely on numerical methods to understand physical behavior.This paper presents a numerical method for chaotic problems.Using fractalfractional differentiation and integral operators in the sense of Newton interpolation polynomial, we investigated the Burke-Shaw system (1) of mathematical equations able to capture chaotic behavior.Solutions are obtained for the fractional-order Burke-Shaw system (1) using a fractional operator with a non-singular kernel.Uniqueness and boundedness for solution are proved through fixed point theory.Due to the high non-linearity of our problem, we used a suitable numerical scheme to solve this system of equations numerically.The presented scheme is applicable to many other systems, see for example, [ 10 , 20 , 24 , 25 ].For similar numerical results, see [41][42][43][44][45][46][47][48][49][50][51].In future work, the existence and uniqueness of solutions reported for general component differential equations will be extended to multidimensional problems.

Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.